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00 . Abstract 

CN ' Deconfinement phase transition due to disappearance of confin- 

ing colorelectric field correlators is described using nonperturbative 
, equation of state. The resulting transition temperature T c (/x) at any 

chemical potential \i is expressed in terms of the change of gluonic con- 
r S^ • densate AG2 and absolute value of Polyakov loop Lt un ^(T c ), known 

from lattice and analytic data, and is in good agreement with lattice 
Ph! data for AG 2 ~ 0.0035 GeV 4 . E.g. T c (0) = 0.27; 0.19; 0.17 GeV for 

rif = 0, 2, 3 respectively. 

> 

X 1 Introduction 



1. Phase transition at nonzero \x and dynamics of quark gluon plasma (qgp) is 
now of great interest because of impressive results of heavy ion experiments, 
see [1] for a recent review and references. The topic calls for a nonperturba- 
tive (NP) treatment of QCD degrees of freedom at nonzero T and /1, which 
is especially important for not very large T, /1. Below we are using the NP 
approach based on Field Correlator Method (FCM) [2], which was applied 
to nonzero T in (31 H] . 

The main advantage of FCM is a natural explanation and treatment of 
dynamics of confinement, as well as the deconfinement transition [31 H], in 
terms of Color Electric (CE) D E {x), Df(x) and Color Magnetic (CM) Gaus- 
sian (quadratic in F" ) correlators D H (x), D^(x). 
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The correlators D E (x) and D H (x) ensure confinement in the planes (4, i) 
and (i,k) respectively, i, k = 1,2,3 so that standard string tension a E = 
a = | J D E (x)d 2 x, and spatial string tension a H = a s = | J D H (x)d 2 x. 
Correlators Df , contain perturbative series, and Df plays an important 
role in that it contributes to the modulus of the Polyakov line; at T > T c one 
has [5J 

4w = \^-(trPexpigj\,dz,)\ = exp (-^^Tj (1) 
with the static QQ potential 

V 1 (r,T)= [ du(l-vT) T^f(v / e+^)- (2) 
Jo Jo 

Deconfinement phase transition in this language is the disappearance of 
D E (x) (and a E ) at T > T c , while L^J nd and L% = (L^ nd ) 9/4 are nonzero 
therej. 

However the disappearance of D E (x) = ■^-{trE i (x)^(x,y)E i (y)) implies 
vanishing of a part of vacuum energy density, 

e vac = ^<(i^) 2 > = - (U ~J nf) G^ D E (0)+D E (0)+D H (0)+D?(0) = ^G 2 

(3) 

and G 2 is the gluon condensate [6]. At T = 0, D E = D H ,D E = £>f and 
for T > both D H , do not change till T « 2T C , while D E disappears at 
T ~> T c [7] in agreement with the deconfinement mechanism suggested in [1] . 

Particle data [7] and analytic study [8] imply that d[ E ' H \x) « 0.2D^-^(x), 
therefore one expects that AG 2 = G 2 (T < T c ) - G 2 (T > T c ) ss |G 2 (T < 
r c ) « |G 2 *, where G 2 * « 0.012 GeV 4 [6] (see [9] for a recent gauge-string 
duality treatment of yielding G 2 * ^ (0.01 ± 0.002) GeV 4 ). 

This AG 2 taken as the change of free energy (pressure) across the phase 
boundary will be our basic element in finding the phase transition curve T c (//) 
below. 

1 The subscript V in t^J^ nd is to distinguish from L^ nd calculated on the lattice with 
singlet free energy F^q(oo,T) replacing Vi(oo,T) in (JlJ. It is clear that L^ F > contains all 
bound states in addition to the ground state (sclfencrgy) Vi(oo,T), hence Fh~. < V\. We 
shall ignore the difference in the first approximation in what follows and write Lf un d- 
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To this end we introduce in the next section the NP EoS of qgp derived 
recently in [JO], and express T c (/x) in terms of AG 2 and Lf und (T c ). 

Taking for the latter the lattice or analytic value, one obtains a set of 
curves T c (/i) for rif = 0,2,3 depending on the only parameter AG2. These 
resulting curves and their end points T c (0), /x c (0) are discussed in conclusion. 

2. In the NP approach to the qgp in [TO] one introduces in the first 
approximation the interaction of single quarks and gluons with the vacuum, 
which is called the Single Line Approximation (SLA), leaving pair and triple, 
etc... correlations to the next steps. As a result one obtains in SLA the 
pressure P^ LA of quarks (and antiquarks) and Pg LA of gluons which are 
expressed through Lf und , namely [TU] . 
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In (BJ, (El) it was assumed that T < y = 1 GeV, where A is the vacuum 

e~\ x \l x , hence powers of L™, see [10] for 



correlation length, e.g. .D^ 
details. 

With few percent accuracy one can replace the sum in (j5j) by the first 
term, n = 1, and this form will be used below for p g i, while for p q this re- 
placement is not valid for large ^, and one can use instead the form equivalent 
to (ID, 



Pq 



TT 2 





(7) 



where v = m q /T and 



z 4 dz 



V^ 2 + v 2 (e^ T + UI ~ a + I) 



(8) 



Eqs. (J7|), (j5J) define p q ,p g i for all T, and m g , which is the current (pole) 
quark mass at the scale of the order of T. 
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Using flU-flS} we can define the pressure Pj in the confined phase, and 
Pu in the deconfined phase, taking into account that vacuum energy density 
in two phases e vac and respectively contributes to the free energy, and 
hence \e vac \, |£*cl to the pressure. One has 

Pi = \^vac\ + Phadron, Pi I = \^ va %\ + Oty + Pq)^. (9) 

Here Phadron is the pressure of the hadron gas at T < T c . From Pi(T c ) = 
Pn(T c ) one obtains T c , neglecting Phadron in the first approximation 

, x /(ll- |n f )AG 2 \ 1/4 , s 



32(Pgi +P q ) 

In (ITOl) enter only two parameters; AG2 and Lf un d{T c ) = exp ^— ^ 
«= ^1(00, T c ) ^^ Q {oo,T c ). 

The latter can be found in 3 different ways: 1) from the direct lattice 
measurements [UJ of Pqq ~ 0.5 GeV; 2) from analytic calculation of Df in 

[8], which yields 14(oo,T < T c ) w 6 " s ^ o)<T « 0.5 GeV with M = (0.8 - 1) 
GeV lowest gluelump mass [12] ; 3) from lattice calculations of -Df at T > T c , 
[H [13], which according to (|2J) yields Vi(oo, T c ) rs 0.4 4- 0.6 GeV. Therefore 
one can fix Vi(oo,T c ) = 0.50(5) GeV (k = 0.25 GeV) and this value is 
independent of rif [IT]. As a result T c (/x) is a function of only AG2 and 
for each value of AG2 one finds a set of curves for n/ = 2,3, ... We choose 
AG 2 » |Gf and in Fig.l the curves computed numerically from ffTUj) for 

= 2, 3 are shown for AG2 ~ 0.00341 GeV 4 and zero quark pole masses. 

The end points T c (0) and /x c (0) can be found analytically e.g. for T c (/x) 
with ~ 5% accuracy one has (expanding ffTOl in ^) 

r c (/i) = T c (0) ( 1 - C-^rV C = 0.0110(3), (11) 



2*(0) 



with 



^°)-2 T() ( 1 + V 1 + ^))' ^ = ( % J • (12) 

For AG 2 = 0.00341 GeV 4 one obtains T c (0) = (0.27; 0.19; 0.17) GeV for 
n/ = 0,2,3 respectively, which agrees well with numerous lattice data, see 
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Figure 1: The phase transition curve T c (n) from Eq. ffTOl) (in GeV) as function 
of quark chemical potential fi (in GeV) for rif = 2 (upper curve) and rif = 3 
(lower curve) and AG 2 = 0.0034 GeV 4 . 



[T4] for reviews. The value of C — 0.011 is inside the scattered set of lattice 
values [T4] . 

Another end point, /x c (0) can be found from the asymptotics of (jHJ), 
$o(a -> oo) = ^ + ^a 2 + ..., which yields (for m g = 0) /i c (0) = k + (48) 1 / 4 T(°) 
and for the same AG 2 as above one gets fi c (0) = (0.63; 0.58) GeV for rif = 2, 3 
One can check, that the derivative in T, vanishes at T = 0. 

3. The phase curve T c (/a) in Fig.l is in reasonable agreement with lattice 
data at least for \i < 0.25 GeV, see [T5] for review and references. Two 
important points are to be discussed here: 1) order of transition and possible 
critical point 2) approximations and assumptions of the present work. 

1). The vacuum transition of our approach is evidently of the first order 
at least in the leading (SLA) approximation used for ffTUl) . and does not 
contain any critical points. This is in agreement with lattice rif — data, 
but the lattice results for rif = 2, 3 depend on masses, discretization and are 
not fully conclusive. The softening of transition for rif > in our approach 
is explained by the increasing role of Phadron{T) for rif > near T c , which 
suppresses the specific heat and makes the curve P(T) more smooth. The 
chiral transition in our approach is caused by the deconfinement, since both 
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(qq) and f n are expressed via D E (x) [T5] and vanish together with it, in 
agreement with lattice data, see e.g. [T4]. The Polyakov loop is a good 
(approximate) order parameter for rif = 0(nf > 0) since at T < T c it is 
expressed via D E (x) and vanishes (strongly decreases) (see Eq. (Q of [5]). 

2) In our derivation of f|T0|) - f|T2|) it was assumed a) that the only important 
part of qgp dynamics is the interaction with the NP vacuum -SLA; b) it is 
assumed that vacuum fields do not depend on T, \x in the phase diagram, 
except at the phase boundary where the shift AG 2 occurs; in particular 
neither AG 2 nor Lf un d(T c ) depend on \i. The latter point is partly supported 
by lattice data [17] . In general this picture of rigid vacuum is based on the 
notion of the dilaton scale of vacuum fields, which can be associated 
with the ++ glueball mass around 1.5 GeV and therefore for all external 
parameters (like ix or T) much less than rrid vacuum fields are fixed. 

Another argument in favor of rigid vacuum is that all dependence on ix and 
rif does not appear in the lowest order of l/iV c expansion, since it comes from 
the quark loops. (Note that nevertheless T c (0) differs strongly for rif = and 
rif = 2,3; even through AG 2 was kept fixed, and this successful prediction 
of T c (0) = 0.27 GeV and 0.19 GeV respectively can be considered as another 
support of our picture). Several things were not taken into account. Quark 
masses are included trivially via Eqs. (161171) and this can be checked vs lattice 
data. Phase transition near /x c (0) can be complicated due to strong qq and 
qqq interaction, which is not taken into account above and will be discussed 
elsewhere (see also [1 Oj ) . therefore the possibility of color superconductivity 
is not commented here. 
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